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Abstract. Although acoustics is one of the disciplines of mechanics, its "geometrization" is still lim- 
ited to a few areas. As shown in the work on nonlinear propagation in Reissner beams, it seems that 
an interpretation of the theories of acoustics through the concepts of differential geometry can help to 
^^ . address the non-linear phenomena in their intrinsic qualities. This results in a field of research aimed 

"^H ■ at establishing and solving dynamic models purged of any artificial nonlinearity by taking advantage of 

^^ , symmetry properties underlying the use of Lie groups. The geometric constructions needed for reduction 

are presented in the context of the "covariant" approach. 

^. 1 Introduction 

^ , The Reissner beam is one of the simplest acoustical system that can be treated in the context of mechanics 

Swith symmetry. A Lie group is a mathematical construction that handle the symmetry but it is also a manifold 
on which a motion can take place. As emphasized by Arnold [1], physical motions of symmetric systems 
governed by the variational principle of least action correspond to geodesic motions on the corresponding 
group G. This paper will try, in a first part, to illustrate this basic concept in the case of the continuous 
^^ ■ group of motion in space. After a literature survey on this subject, an extension from geodesies to auto- 

^T \ parallel submanifolds is proposed in the second part and naturally leads to the geometric covariant approach 

Q^ ■ available to study evolution problems for fields defined by a variational principle. 

p 

"^ ! 2 Nonlinear model for Reissner Beam 

o_ 

^^ ! 2.1 Reissner kinematics 



A beam of length L, with cross-sectional area A and mass per unit volume p is considered. Following the 
Reissner kinematics, each section of the beam is supposed to be a rigid body. The beam configuration can 
/\^ . be described by a position r(^, t) and a rotation R(^, t) of each section. The coordinate s corresponds to the 

position of the section in the reference configuration i7o (see figure [TJ. 



2.2 Lie group configuration space 

Any material point M of the beam which is located at x(^, 0) = r(s, 0) + wq = ^Ei + wq in the reference 
configuration (t = 0) have a new position (at time t) x(s, t) = r(s, t) + R(5', Owq. In other words, the current 
configuration of the beam Zt is completely described by a map 

x(M)|^|R(jOr(.,0||wo|^ R.50(3), r.R3, (1) 




Fig. 1: Reference and current configuration of a beam. Each section, located at position s in the reference 
configuration i7o, is parametrized by a translation r(^, t) and a rotation R(^, t) e SO3 in the current 

configuration i7^ 



where the matrix H(s, t) is an element of the Lie group SE{?>) = S0(3) x R-^, where S0(3) is the group 
of all 3 X 3 orthogonal matrices with determinant 1 (rotation in R^). As a consequence, to any motion 
of the beam a function H(^, t) of the (scalar) independent variables s and t can be associated. Given some 
boundary conditions, among all such motions, only a few correspond to physical ones. What are the physical 
constraints that such motions are subjected to? 

In order to formulate those constraints the definition of the Lie algebra is helpful. To every Lie group 
G, we can associate a Lie algebra g, whose underlying vector space is the tangent space of G at the identity 
element, which completely captures the local structure of the group. Concretely, the tangent vectors, ^^H 
and dtH, to the group S E(3) at the point H, are lifted to the tangent space at the identity e of the group. 
The definition in general is somewhat technical , but in the case of matrix groups this process is simply a 
multiplication by the inverse matrix H"^ This operation gives rise to definition of two left invariant vector 
fields in g = 6e(3) 






(2) 
(3) 



which describe the deformations and the velocities of the beam. Assuming a linear stress-strain relation, 
those definitions allow to define a reduced Lagrangian function by the diff'erence of kinetic and potential 



^ In the literature, one can find the expression dL^- 
an element of G defined by 



(g) where dL stands for the differential of the left translation L by 



L,:G 



G 

hog. 



energy K^^, €c) = Ec - Ep, with 



Ec(Xc) = f ^x'JXcds, (4) 

Ep(6c)= I ^(ec-eofC(6c-eo)ds, (5) 



2^c-^c- 



^^^=X^^^^ 



where J and C are matrix of inertia and Hooke tensor respectively and 60 = H ^(s, 0)^^H(^, 0) correspond 
to the deformation of the initial configuration. 

2.3 Equations of motion 

Applying the Hamilton principle to the left invariant Lagrangian / leads to the Euler-Poincare equation 

dtTTc - cid^^TTc = ds{(Tc - ctq) - adl{(Tc - o-q), (6) 

where Tic = JXc ^^^ ^c = ^6^, (see for example fT|, f4l or f5l for details). In order to obtain a well-posed 
problem, the compatibility condition, obtained by diff'erentiating Q and (O 

dsXc - ^t^c = ad^^Sc, (7) 

must be added to the equation of motion. It should be noted that the operators ad and ad"" in eq. Q 

ad*^^ v)("^' p) = (m X a> + p X V, p X a>) (8) 

cid^oji,Yi)(co2, V2) = (coi X a>2, a>i X V2 - a>2 X vi), (9) 

depend only on the group S E(3) and not on the choice of the particular "metric" L that has been chosen to 
described the physical problem [6|. 

Equations © and d?]) are written in material (or left invariant) form (c subscript). Spatial (or right invari- 
ant ) form exist also. In this case, spatial variables (s subscript) are introduced by 

es(s,t) = dsH(s,m-\s,t) (10) 

X\(sj) = dtH(s,m-\sj) (11) 

and Q leads to the conservation law [18 1 

dtTls = dsiCTs - CTo) (12) 

where tTs = Ad^.i tt^ and o"^ = Adg_iO"c. The Ad* map for S E(3) is 

Adg 1 (m, p) = (Rm + r x Rp, Rp). (13) 

Compatibility condition ^ becomes 

dsXs - ^t^s = ade^Xs- (14) 

Equations Q and ^ (or alternatively ([12]) and ([14])) provide the exact non linear Reissner beam model and 
can be used to handle the behavior of the beam if the large displacements are taking into account. 

Notations and assumptions vary so much in the literature, it is often difficult to recognize this model (see 
for example |7| for a formulation using quaternions). However, this generic statement is used to classify 



publications according to three axes. In the first one, the geometrically exact beam model is the basis for 
numerical formulations. Starting with the work of Simo [2 1, special attention is focused on energy and mo- 
mentum conserving algorithms [8|, [9|. Numerical solutions for planar motion are also investigated in [10|. 
Even, in some special sub-cases (namely where the longitudinal variables do not appear) the non-linear beam 
model gives rise to linear equations which can be solved by analytical methods ifTTIl . 

Secondly, much of the literature is also devoted to the so-called Kirchhoff^'s rod model. In this case, shear 
strain is not taken into account along a thin rod (i.e., its cross-section radius is much smaller than its length 
and its curvature at all points). In this approximation cross-sections are perpendicular to the central axis of 
the filament and the rotation matrix can be given in the Frenet-Serret frame, (see [12] , [14], [15], for 
example). In that context an interesting geometric correspondence between Kirchhoff'rod and Lagrange top 
can be made |13|. 

Finally, if only rigid motion is investigated, ( i.e. if the spatial dependence in ^ is canceled: ds = 0) 
the so-called underwater vehicleEl model is obtained. In absence of exterior force and torque, the equation of 
motion for a rigid body in an ideal fluid become more simply [|T6ll . ifTTl 

^ |m = mX6L>-hpxv 

dtTTc = ad^ Tie, that is < (15) 

In this simpler form, a geometric interpretation is easier. The solution of the equation of motion mentioned 
above, if it exists, should be interpreted as a geodesic of the group S E(3) endowed with a non-canonical left 
invariant metric J. To accomplish the correspondence between the Euler-Poinc are's equation and geodesic 
equation the historical definition of the covariant derivative is exposed in the next section. 

3 Geometric interpretation 

3.1 Geodesies on curved spaces 

A trajectory of a particle of mass m which is moving on a manifolQj M can be thought as a curve a(t) 
on M and v(0 = cx(t) is the speed of the particle. According to the Newton's second Law of motion, its 
acceleration (the variation of its velocity) is proportional to the net force acting upon it 2 F = m^. The 
expression of this variation, v(^ -h dt) - v(0, shows that the velocities are evaluated at two diff'erent points 
of the curve: a(t -\- dt) and a(t) which are, a priori, incommensurable quantities. So, one of the two vectors 
needs to be parallel transport as it is illustrated, for flat manifolds, in figure (|2]). For curved manifolds the 




\(t + dt) - \(t) 



Fig. 2: For flat manifolds, a trivial 

parallel transport is used to compute 

the acceleration. 



operation is not so easy and its historical construction is related by M.P. do Carmo in |19| for surfaces of 
R^ (see figure [S]). Technically, this historical construction gives rise to the concept of the covariant derivative 

^ underwater vehicle in the case that the center of buoyancy and the center of gravity are coincident 
^ a surface for short 




d<t>„(W> 
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Fig. 3 Parallel transport along a curve: Let a(t) be 
a curve on a surface S and consider the envelope 
of the family of tangent planes of S along a (see 
figureO. Assuming that a(t) is nowhere tangent to 
an asymptotic direction, this envelope is a regular 
surface i7 which is tangent to S along a. Thus, the 
parallel transport along a of any vector we Tp(S), 
p e S ,isthQ same whether we consider it relative 
to S or to i7. Furthermore, i7 is a developable sur- 
face; hence can be mapped by an isometry into 
a plane P (without stretching or tearing). Parallel 
transport of a vector w is then obtained using usual 
parallel transport in the plane along (p(a) and pull 
it back to i7 (by d(fi~^). 



^ = VyW of a vector field w along a. The parametrized curves a : I ^ R^ of a. plane along which the field 
of their tangent vector v(0 is parallel are precisely the straight lines of that plane. The curves that satisfy an 

analogous condition , i.e. 

D\ 

— = V,v = 0, (16) 

at 

for a surface are called geodesies. Intuitively, the acceleration as seen from the surface vanishes : in absence 

of net force, the particle goes neither left nor right, but straight ahead. 

The kinetic energy dU define a left invariant Riemannian metric on S E(3), and then define also a sym- 
metric connection V which is compatible with this metric (Levi-Civita connection). It can be shown that 
geodesic equation ([T6b for this particular connection coincide with Euler-Poincare equation of motion ([TSll 
when S E(3) is endowed with the kinetic metric (|4]). 

Now, this equation deals with motion of rigid body described by a single scalar variable t. So what is the 
geometric interpretation of the equations of motion ^ and d?]) for which two variables s and t are involved. 
In other word, can we extend a geodesic, which is a 1 -dimensional manifold, to 2-dimensional geodesic ? 



3.2 Auto-Parallel submanifolds, covariant point of view 

A geodesic curve on a surface 5" is a 1 -dimensional submanifold of S for which the parallel transport 
of its initial velocity stay in its own tangent space. In that sense, a geodesic is an auto-parallel curve. If 
now, geodesies are seen as auto-parallel curves on surface, a definition of an n-dimensional auto-parallel 
submanifolds can be made. 

A submanifold M is auto-parallel in S if the parallel translation of any tangent vector of M along any 
curve in M stays in its own tangent space T(M). Note that a parallel translation of a vector w e T(M) 



certainly belongs to T(S) but not necessarily to T(M). In other words, M is auto-parallel in S with respect 
to the connection V of 5" if VxY belongs to T(M), VX, Y g T(M). A correspondence between auto-parallel 




Fig. 4: Symbolic representation of a parametrized surface i7 immersed into the group G = S E(3) 



surface and solutions to equations (O and d?]) is still to be demonstrated. In this case, any motion of the beam 
must be seen as a map from [0, L] x R c R^ to S E(3) given by 



(s, t) -^ His, 0, 



(17) 



rather then a curve t -^ H(., t) in the infinite dimensional configuration space 9^([0; L];SE(3)) of functions 
from [0; L] to S E(3). In this perspective, solving a physical variational problem is therefore transposed to the 
problem of finding an auto-parallel immersed surface as it is illustrated symbolically in figure (|4]). This pro- 
cess illustrates the covariant (as opposed to dynamical) formulation of a variational problem (see [1211 . [|22l ). 
More precisely, the map SVJ\ should be interpreted as a (local) section 6(x) = (x, H(x)) of the principal 
fiber bundle P ^ X with structure group G = S E(3) 

n:P = Xx SE(3) -^ X, 7r(x, H) := x 

over the spacetime X = [0; L] x R, (s, t) = x. The Lagrangian is then defined in the phase bundle L : 
J^P -^ R, where J^P denotes the first jet bundle of the bundle P. If L is invariant under the action of G, 
the variational principle drops to the quotient space (J^P)/G. This quotient is an affine bundle on X which 
can be identified to the bundle of connections C ^ X. It induces a reduced lagrangian I \ C ^ X from L 
and a reduced section s g FiC) from s g r{P). In that context, the multidimensional generalization of the 
equations of motion (|6]), compatibility ^ and conservation law ([T2l) are formulated by M. Castrillon Lopez 
in [231 . The equation of motion (O now yields 



div— + ad^ — 
61 '61 







(18) 



where div stands for the divergence operator defined by the volume form v (here, v = ds Adt). The compat- 
ibility condition ©gives rise to the flatness of s (integrability condition) 



curv(6) = ds-\- [s, s] = 0. 



(19) 



Finally, introducing the Cartan-Poincarqj n-form, Sl, the symmetries of a variational problem produces 
conservation laws by means of the Noether's Theorem 

j(/i5)*J" = 0. (20) 

The form J^ induces a conserved quantity since its differential vanishes along the critical section s. It can be 
understood as a current form (like in electromagnetism). In that sense, this formulation is more appropriated 
to describe a conservation law then the partial derivative balance law ([T2l) . But for a non-specialist audience, 
the definition of this form 

is quite obscure (in particular the relationship between the Cartan-Poincare form and the conserved quantity). 
Here j stands for the interior product and ^* is related to the infinitesimal vector field generated by the 
symmetr}0. 

4 Conclusion 

A geometrical approach of the dynamic of a Reissner beam has been studied in this article in order to 
take into account non linear effects due to large displacements. There are basically two different geomet- 
ric approaches available to study evolution problems for fields defined by a variational principle. The first 
approach, called the "dynamical" approach, uses, as its main ingredient, the infinite dimensional manifold 
as configuration space (TQ). The reduction techniques developed in the dynamical framework have been 
studied thoroughly in the literature (see for example [24] and the references therein cited), but it presents the 
difficulty to handle geodesic curves in an infinite dimensional function space. 

As an alternative, the covariant formulation allows to consider a finite dimensional configuration space 
(the dimension of the symmetry group itself). Although its roots go back to De Bonder |25|, Weyl |^6l, 
Caratheodory [27 1, after J. M. Souriau in the seventies [28 1, the classical field theory has been only well un- 
derstood in the late 20th century (see for example [29] for an extension from symplectic to multisymplectic 
form). It is therefore not surprising that, in this covariant or jet formulation setting, the geometric construc- 
tions needed for reduction have been presented even more recently. In that circumstances in the literature, 
it is also not easy to understand how the multisymplectic form can be obtained from the differential of the 
Cartan-Poincare n-form, which is crucial to give rise to an Hamiltonian framework (Lie-Poisson Schouten- 
Nijenhuis (SN) brackets [30]). An understandable theory, that can unify all the results obtained "ad hoc", 
case by case, is still missing to our knowledge. 
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